Thermalization and its mechanism for generic isolated quantum systems* 
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Time dynamics of isolated many-body quantum systems has long been an elusive subject. Very recently, 
however, meaningful experimental studies of the problem have finally become possible [HQ], stimulating the- 
oretical interest as well H,0>U,@,0]. Progress in this field is perhaps most urgently needed in the foundations 
of quantum statistical mechanics. This is so because in generic isolated systems, one expects 0| @] nonequi- 
librium dynamics on its own to result in thermalization: a relaxation to states where the values of macroscopic 
quantities are stationary, universal with respect to widely differing initial conditions, and predictable through the 
time-tested recipe of statistical mechanics. However, it is not obvious what feature of many-body quantum me- 
chanics makes quantum thermalization possible, in a sense analogous to that in which dynamical chaos makes 
classical thermalization possible 1 10]. For example, dynamical chaos itself cannot occur in an isolated quantum 
system, where time evolution is linear and the spectrum is discrete 1 11]. Underscoring that new rules could 
apply in this case, some recent studies even suggested that statistical mechanics may give wrong predictions 
for the outcomes of relaxation in such systems |4j,|5|]. Here we demonstrate that an isolated generic quantum 
many-body system does in fact relax to a state well-described by the standard statistical mechanical prescription. 
Moreover, we show that time evolution itself plays a merely auxiliary role in relaxation and that thermalization 
happens instead at the level of individual eigenstates, as first proposed by J. M. Deutsch 1 12] and M. Srednicki 
11311 - A striking consequence of this eigenstate thermalization scenario, confirmed below for our system, is that 
the knowledge of a single many-body eigenstate suffices to compute thermal averages — any eigenstate in the 
microcanonical energy window will do, as they all give the same result. 



If we pierce an inflated balloon inside a vacuum chamber, 
very soon we find the released air uniformly filling the en- 
closure and the air molecules attaining the Maxwell velocity 
distribution whose width depends only on their total number 
and energy. Different balloon shapes, placements, or piercing 
points all lead to the same spatial and velocity distributions. 
Classical physics explains this thermodynamical universality 
as follows [10]: almost all particle trajectories quickly begin 
looking alike, even if their initial points are very different, be- 
cause nonlinear equations drive them to explore ergodically 
the constant-energy manifold, covering it uniformly with re- 
spect to precisely the microcanonical measure. However, if 
the system possesses further conserved quantities /Mncf/onaWy 
independent from the Hamiltonian and each other, then time 
evolution is confined to a highly restricted hypersurface of the 
energy manifold. Hence, microcanonical predictions fail and 
the system does not thermalize. 

On the other hand, in isolated quantum systems not only is 
dynamical chaos absent due to the linearity of time evolution 
and the discreteness of spectra [11], but it is also not clear un- 
der what conditions conserved quantities provide independent 
constraints on relaxation dynamics. To begin with, any opera- 
tor commuting with a generic and thus nondegenerate Hamil- 
tonian is functionally dependent on it |3, implying that the 
conservation of energy is the only independent constrain. On 
the other hand, even when operators are functionally depen- 
dent, their expectation values — considered as functionals of 
states — generally are not: for example, two states may have 
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the same mean energies but different mean square-energies. 
For nondegenerate Hamiltonians a maximal set of constants 
of motion with functionally independent expectation values is 
as large as the dimension of the Hilbert space; examples in- 
clude the projectors P a — \^ a ) | to the energy eigenstates 
[14] and the integer powers of the Hamiltonian [3[. 

The current numerical and analytic evidence from inte- 
grable systems suggests that there exists a minimal set of in- 
dependent constraints whose size is much smaller than the di- 
mension of the Hilbert space but may still be much greater 
than one. In our previous work (3[[ (with V. Yurovsky) we 
showed that an integrable isolated one-dimensional system of 
lattice hard-core bosons relaxes to an equilibrium character- 
ized not by the usual but by a generalized Gibbs ensemble. In- 
stead of just the energy, the Gibbs exponent contained a linear 
combination of conserved quantities — the occupations of the 
eigenstates of the corresponding Jordan-Wigner fermions — 
whose number was still only a tiny fraction of the dimension 
of the Hilbert space. Yet this ensemble works, while the usual 
one does not, for a wide variety of initial conditions [15] as 
well as for a fermionic system 01611 : it also explains a re- 
cent experimental result, the absence of thermalization in the 
Tonks-Girardeau gas [ 1 ] . Thus, while at least some constraints 
beyond the conservation of energy must be kept, it turns out 
one needs only a relatively limited number of additional con- 
served quantities with functionally independent expectation 
values; adding still further ones is redundant. 

Since it is not clear which sets of conserved quantities — 
and some are always present — constrain relaxation and which 
do not, it becomes even more urgent to determine whether 
isolated generic quantum systems relax to the usual thermal 
state. The theoretical attention to this question has in fact 
been increasing recently, because of the high levels of iso- 
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FIG. 1: Relaxation dynamics, a, Two-dimensional lattice on which 
five hard-core bosons propagate in time. The bosons are initially 
prepared in the ground state of the sub-lattice in the lower-right cor- 
ner and released through the indicated link, b, The corresponding 
relaxation dynamics of the marginal momentum distribution center 
[n(k x = 0)] compared with the predictions of the three ensembles. 
In the microcanonical case, we averaged over all eigenstates whose 
energies lie within a narrow window (see Supplementary Discussion) 
[E - AE,E + AE], where E = (ip(0)\H\ip(0)) = -5.06J, 
AE = 0.1J, and J is the hopping parameter. The canonical en- 
semble temperature is k^T = 1.87J, where /cb is the Boltzmann 
constant, so that the ensemble prediction for the energy is Eo. c, Full 
momentum distribution function in the initial state, after relaxation, 
and in the different ensembles. Here d is the lattice constant and 
L x = 5 the lattice width. 



lation HI HI 03] and control |18, 19] possible in experiments 
with ultracold quantum gases. However, despite numerous 
studies of specific models there is not yet consensus on how 
or even if relaxation to the usual thermal values occurs for 
nonintegrable systems Q. Common wisdom says that it does 
iS IH, but some recent numerical results suggest otherwise, 
either under certain conditions |4J or in general |5y. 

In order to study relaxation of an isolated quantum sys- 
tems, we considered the time evolution of five hard-core 
bosons with additional weak nearest-neighbor repulsions, on 
a 21 -site two-dimensional lattice, initially confined to a por- 
tion of the lattice and prepared in their ground state there. 
Figure QJ shows the exact geometry (see also Supplemen- 



tary Discussion); the relaxation dynamics begins when the 
confinement is lifted. Expanding the initial state wave- 
function in the eigenstate basis of the final Hamiltonian H 
as IV'(O)) = Yla C"a\^a), the many-body wavefunction 
evolves as \tp(t)) = e~ ij& * |-0(O)) = C a e-* E <**|# a ), 
where the E a 's are the eigenstate energies. Thus obtaining 
numerically-exact results for all times required the full di- 
agonalization of the 20,349-dimensional Hamiltonian. The 
quantum-mechanical mean of any observable A evolves as 

(A(t)) ee (V(t)|%(t)> =J2 C aC^- E ^A a0 , (1) 



with A aj 3 
then 



(ty a >\A\^ a ). The long-time average of (A(t)) is 



(A) — ^ ' A aa . 



(2) 



Note that if the system relaxes at all, it must be to this value. 
We find it convenient to think of Eq. (O as stating the pre- 
diction of a "diagonal ensemble," \C a \ 2 corresponding to the 
weight \^ a ) has in the ensemble. In fact, this ensemble is 
precisely the generalized Gibbs ensemble introduced in Ref. 
Q3f] if as integrals of motion one takes all the projection opera- 
tors P a — \^ a ) (*&a\- Using these as constraints on relaxation 
dynamics, the theory gives p c = cxp Y^a=i ^Paj , with 

A Q = — lndCcj 2 ), and D the dimension of the Hilbert space. 
(Notice, however, that for the integrable system treated in Ref. 
||3|], the generalized Gibbs ensemble was defined using a dif- 
ferent, minimal set of independent integrals of motion, whose 
number was equal to the number of lattice sites N -C D.) 

Now if the quantum-mechanical mean of an observable be- 
haves thermally it should settle to the prediction of an appro- 
priate statistical-mechanical ensemble. For our numerical ex- 
periments we chose to monitor the marginal momentum dis- 
tribution along the horizontal axis n(k x ) and its central com- 
ponent n(k x = 0) (see Supplementary Discussion). Figures 
QJ> andQJ demonstrate that both relax to their microcanonical 
predictions. The diagonal ensemble predictions are indistinct 
from these, but the canonical ones, although quite close, are 
not. This is an indication of the relevance of finite size effects, 
which may be the origin of some of the apparent deviations 
from thermodynamics seen in the recent numerical studies of 
Refs. H and M- 

The statement that the diagonal and microcanonical ensem- 
bles give the same predictions for the relaxed value of A reads 



^ I C a | A aa — (^4-)microcan. 



def. 



1 



^ A aa j (3) 



\E -E a \<AE 

where Eq is the mean energy of the initial state, AE is the 
half-width of an appropriately chosen (see Supplementary 
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FIG. 2: Thermalization in classical vs quantum mechanics, a, In 

classical mechanics, time evolution constructs the thermal state from 
an initial state that generally bears no resemblance to the former, 
b, In quantum mechanics, according to the eigenstate thermaliza- 
tion hypothesis, every eigenstate of the Hamiltonian always implic- 
itly contains a thermal state. The coherence between the eigenstates 
initially hides it, but time dynamics reveals it through dephasing. 



Discussion) energy window centered at Eq, and the normal- 
ization Me .ae is the number of energy eigenstates with en- 
ergies in the window [Eq — AE, E + AE}. Thermodynam- 
ical universality is evident in this equality: while the left hand 
side depends on the details of the initial conditions through 
the set of coefficients C a , the right hand side depends only on 
the total energy, which is the same for many different initial 
conditions. Three mechanisms suggest themselves as possi- 
ble explanations of this universality (assuming the initial state 
is sufficiently narrow in energy, as is normally the case — see 
Supplementary Discussion): 

(i) Even for eigenstates close in energy, there are large 
eigenstate-to-eigenstate fluctuations of both the eigenstate ex- 
pectation values A aa and of the eigenstate occupation num- 
bers |C Q | 2 . However, for physically interesting initial condi- 
tions, the fluctuations in the two quantities are uncorrelated. 
A given initial state then performs an unbiased sampling of 
the distribution of the eigenstate expectation values A aa , re- 
sulting in Eq. (O. 

(ii) For physically interesting initial conditions, the eigen- 
state occupation numbers \C a | 2 practically do not fluctuate at 
all between eigenstates that are close in energy. Again, Eq. 
(0 immediately follows. 

(iii) The eigenstate expectation values A aa practically do 
not fluctuate at all between eigenstates that are close in energy. 
In that case Eq. (O holds for literally all initial states narrow 
in energy. 

J. M. Deutsch and M. Srednicki have independently pro- 
posed the scenario (iii), dubbed the 

Eigenstate thermalization hypothesis (ETH) 
[Deutschlflll (1991), Srednicki IB 1 (1994)1. 
The expectation value (>I' ,| J 4| , I' Q ,) of a few-body 



observable A in an eigenstate of the Hamiltonian \^> a ), 
with energy E a , of a large interacting many -body 
system equals the thermal (microcanonical in our case) 
average (A) m i aoclln (E a ) of A at the mean energy E a : 



(*a|A|* a ) = (A) 



microcan. 



(E a ). 



(4) 



The ETH suggests that classical and quantum thermal states 
have very different natures, as depicted in Fig. [2] While at 
present there are no general theoretical arguments supporting 
the ETH, some results do exist for restricted classes of sys- 
tems. To begin with, Deutsch 111 211 showed that the ETH holds 
in the case of an integrable Hamiltonian weakly perturbed 
by a single matrix taken from a random Gaussian ensemble. 
Next, nuclear shell model calculations have shown that indi- 
vidual wavefunctions reproduce thermodynamic predictions 
ll20ll . Then there are rigorous proofs that some particular 
quantum systems, whose classical counterparts are chaotic, 
satisfy the ETH in the semiclassical limit lull 



23, 24]. 



More generally, for low density billiards in the semiclassical 
regime, the ETH follows from Berry's conjecture lfl3i l25ll . 
which in turn is believed to hold in semiclassical classically- 
chaotic systems lEill . Finally, at the other end of the chaos- 
integrability continuum, in systems solvable by Bethe ansatz, 
observables are smooth functions of the integrals of motion. 
This allows the construction of single energy eigenstates that 
reproduce thermal predictions 12711 . 

In Figs. [3^-c we demonstrate that the ETH is in fact the 
mechanism responsible for thermal behavior in our noninte- 
grable system. Fig. [3J: additionally shows that scenario (ii) 
mentioned above plays no role, because the fluctuations in the 
eigenstate occupation numbers \C a \ 2 are large. Thermal be- 
havior also requires that both the diagonal and the chosen ther- 
mal ensemble have sufficiently narrow energy distributions 
p(E) [ = probability distribution x the density of states], so 
that in the energy region where the energy distributions p(E) 
are appreciable, the derivative of the curve eigenstate expecta- 
tion value A aa vs the energy (here n(k x — 0) vs the energy) 
does not change much; see Supplementary Discussion. As 
shown in Fig.|3j), this holds for the microcanonical and diago- 
nal ensembles but not for the canonical ensemble, explaining 
the failure of the latter to describe the relaxation in Fig. Q] 
Note that the fluctuations of the eigenstate occupation num- 
bers \C a \ 2 in Fig. |3]5 are lowered by the averaging involved 
in the computation of the density of states (compare with Fig. 

St). 

To strengthen the case for the ETH, we tested another ob- 
servable. We chose it with the following consideration in 
mind: in our system interactions are local in space, and mo- 
mentum distribution is a global, approximately spatially addi- 
tive property. Thus one might wonder if the ETH for momen- 
tum distribution arises through some spatial averaging mech- 
anism (we thank the anonymous referee 2 for bringing this 
point to our attention). It does not: for our final test of the ETH 
we chose an observable that is manifestly local in space, the 
expectation value of the occupation of the central site of the 



4 



lattice. We again found that the ETH holds true (3% relative 
standard deviation of eigenstate-to-eigenstate fluctuations). 

On the other hand, Figs. [3}i-f show how the ETH fails 
for an isolated one-dimensional integrable system. The lat- 
ter consists of five hard-core bosons initially prepared in their 
ground state in an 8-site chain, one of the ends of which we 

A noniniegrable 



then link to one of the ends of an adjoining (empty) 13-site 
chain to trigger relaxation dynamics. As Fig. [3^ shows, n(k x ) 
as a function of energy is a broad cloud of points, meaning 
that the ETH is not valid; Fig. [3j shows that scenario (ii) does 
not hold either. 
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FIG. 3: Eigenstate thermalization hypothesis, a, In our nonintegrable system, the momentum distribution n(k x ) for two typical eigenstates 
with energies close to Eq is identical to the microcanonical result, in accordance with the ETH. b, Upper panel: n(k x — 0) eigenstate 
expectation values as a function of the eigenstate energy resemble a smooth curve. Lower panel: the energy distribution p(E) of the three 
ensembles considered in this work, c, Detailed view of n(k x — 0) (left labels) and \C a \ 2 (right labels) for 20 eigenstates around Eq. d, In the 
integrable system, n(k x ) for two eigenstates with energies close to Eo and for the microcanonical and diagonal ensembles are very different 
from each other, i.e., the ETH fails, e, Upper panel: n(k x — 0) eigenstate expectation value considered as a function of the eigenstate energy 
gives a thick cloud of points rather than resembling a smooth curve. Lower panel: energy distributions in the integrable system are similar to 
the nonintegrable ones depicted in b. f, Correlation between n(k x = 0) and |C a | 2 for 20 eigenstates around Eq. It explains why in d the 
microcanonical prediction for n(k x — 0) is larger than the diagonal one. 



Nevertheless, one may still wonder if in this case scenario 
(i) might hold — if the averages over the diagonal and the 
microcanonical energy distributions shown in Fig. [3^ might 
agree. Figure [3}i shows that this does not happen. This is so 
because, as shown in Fig. [3f, the values of n(k x = 0) for 
the most-occupied states in the diagonal ensemble (the largest 
values of eigenstate occupation numbers |Co,| 2 ) are always 
smaller than the microcanonical prediction, and those of the 
least-occupied states, always larger. Hence, the usual thermal 
predictions fail because the correlations between the values 
of n(k x = 0) and \C a \ 2 preclude unbiased sampling of the 
latter by the former. These correlations have their origin in 
the nontrivial integrals of motion that make the system inte- 



grable and that enter the generalized Gibbs ensemble, which 
was introduced in Ref. yfl as appropriate for formulating sta- 
tistical mechanics of isolated integrable systems. In the non- 
integrable case shown in Fig. [3J;, n(k x = 0) is so narrowly 
distributed that it does not matter whether or not it is corre- 
lated with | C a | 2 (we have in fact seen no correlations in the 
nonintegrable case). 

The thermalization mechanism outlined thus far explains 
why long-time averages converge to their thermal predictions. 
A striking aspect of Fig. QJ), however, is that the time fluc- 
tuations are so small that after relaxation the thermal predic- 
tion works well at every instant of time. Looking at Eq. (fl}, 
one might think this is so because the contribution of the off- 
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FIG. 4: Temporal vs quantum fluctuations, a, Matrix elements of 
the observable of interest, n(k x — 0), as a function of state indices; 
the eigenstates of the Hamiltonian are indexed in the order of dimin- 
ishing overlap with the initial state. The dominance of the diagonal 
matrix elements is apparent, b, The same time evolution as in Fig. [Tp 
with the error bars showing the quantum fluctuations n(k x — 0) ± A 
with A = [{n 2 (k x = 0)) - {n(k x = 0)) 2 ] 1/2 , which are clearly 
much larger than the temporal fluctuations of n(k x = 0). 



diagonal terms gets attenuated by temporal dephasing, which 
results from the generic incommensurability of the frequen- 
cies of the oscillating exponentials. However, this attenuation 
only scales as the root of the number of dephasing terms, and 
is exactly compensated by their larger number: if the number 
of eigenstates that have a significant overlap with the initial 
state is -/V states , then typical C a ~ W^Wes, and the sum 
over off-diagonal terms in Eq. (Q} finally does not scale down 
with AWes: 



E 

a, (3 



i(E a -Ef))t 



N s , 



-A 



af3 



states ^typical 



iVst 



a/3 



A 



typical 

a/3 



(5) 



Hence, were the magnitude of the diagonal and off-diagonal 
terms comparable, their contributions would be comparable 
as well, and time fluctuations of the average would be of the 
order of the average. However, this is not the case and thus 



^typical ^typical 

ajt/3 



(6) 



Figure confirms this inequality for the matrix elements of 
the momentum distribution in our system. We should mention 
that there is an a priori argument — admittedly in part depen- 
dent on certain hypotheses about chaos in quantum billiards — 
in support of this inequality for the case when the mean value 
of A in an energy eigenstate is comparable to the quantum 
fluctuation of A in that state [28]. 

On the other hand, the thermalization we see appears to 
be working a bit too well: in a system as small as ours, one 
would expect measurement-to-measurement fluctuations to be 
much larger than what Fig. QJ) suggests. Indeed, as we show 
in Figure^, the fluctuations that one would actually measure 
would be dominated by the quantum fluctuations of the time- 
dependent state. The rather large size of the quantum fluctu- 



ations relative to the thermal mean value is of course particu- 
lar to small systems; however, the dominance of the quantum 
fluctuations over the temporal fluctuations of quantum expec- 
tation values is not and is actually expected for generic sys- 
tems in the thermodynamic limit II29I1 . 

We have demonstrated that, in contrast to the integrable 
case, the nonequilibrium dynamics of a generic isolated quan- 
tum system does lead to standard thermalization. We verified 
that this happens through the eigenstate thermalization mech- 
anism, a scenario J. M. Deutsch [12] demonstrated for the case 
of an integrable quantum Hamiltonian weakly perturbed by a 
single matrix taken from a random Gaussian ensemble and M. 
Srednicki IU3I1 compellingly defended for the case of rarefied 
semiclassical quantum billiards, and which both authors con- 
jectured to be valid in general. Our results, when combined 
with the others we mentioned (12. [13. |23. EH E3. E3. EH 
26l 12711 . constitute strong evidence that eigenstate thermaliza- 
tion indeed generally underlies thermal relaxation in isolated 
quantum systems. Therefore, to understand the existence of 
universal thermal time-asymptotic states, one should study 
operator expectation values in individual eigenstates. This is 
a problem that is linear, time-independent, and conceptually 
far simpler than any arising in the research — currently domi- 
nating the field — on the nonlinear dynamics of semiclassical 
systems. Among the fundamental open problems of statistical 
mechanics that could benefit from the linear time-independent 
perspective are the nature of irreversibility, the existence of a 
KAM-like threshold [30] in quantum systems, and the role of 
conserved quantities in the approach to equilibrium. Finally, 
having a clear conceptual picture for the origins of thermaliza- 
tion may make it possible to engineer new, "unthermalizable" 
states of matter [12], with further applications in quantum in- 
formation and precision measurement. 

We thank A. C. Cassidy, K. Jacobs, A. P. Young, and E. 
J. Heller for helpful comments. We acknowledge financial 
support from the National Science Foundation and the Office 
of Naval Research. We are grateful to the USC HPCC center 
where our numerical computations were performed. 



SUPPLEMENTARY DISCUSSION 



1. The Hamiltonian and the numerical calculations. 



In a system of units where h — 1 the Hamiltonian reads 



(i,3) 



b\bj + h.c 



U fiifij 

(i,3) 



(7) 



where indicates that the sums run over all nearest- 

neighbor pairs of sites, J is the hopping parameter, and U 
the nearest-neighbor repulsion parameter that we always set 
to 0.1J. The hard-core boson creation (Sj) and annihilation 
(bj) operators commute on different sites, [b i , bj] = [b i , bj] = 

[SJ, bj] = for all i and j ^ i, while the hard-core condi- 
tion imposes the canonical anticommutation relations on the 
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same site, {b i , b\} — 1, and = (hi) 2 = for all i. Here 
hi = b\b i is the density operator. 

An exact study of the nonequilibrium dynamics for all time 
scales requires a full diagonalization of the many -body Hamil- 
tonian (0. We are able to fully diagonalize — essentially to 
machine precision — matrices of dimension D ~ 20, 000, and 
so we consider N = 5 hard-core bosons on a 5 x 5 lattice 
with four sites missing (D = 20, 349); see Fig. [5] All the 
eigenstates of the Hamiltonian are used for the time evolution 

\i>(t)) = exp \-iHt] |^(0)) = C « ex P [-iEJ\\*a) , 

a 

where \ip(t)) is the time-evolving state, \ip(0)) the initial state, 
\^ a ) the eigenstates of the Hamiltonian with the energies E a , 
and C a — ($ Q: |'0(O)). Our initial state is the ground state of 
the five bosons when they are confined to the lower part of the 
lattice (the colored part in Fig. [5] The time evolution begins 
with the opening of the link shown in Fig. [5] which allows the 
bosons to expand over the whole lattice. The position of the 
missing sites was chosen so that we only open a single link 
to start the relaxation dynamics. The motivation for this will 
become apparent in the last paragraph below. 



P Initial 
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FIG. 5: The lattice for the dynamics. Two-dimensional lattice on 
which the particles propagate in time. The initial state is the ground 
state of 5 hard-core bosons confined to the sub-lattice in the lower 
right-hand corner, and the time evolution starts after the opening of 
the link indicated by the door symbol. 



within a narrow window [Eq — AE,Eq + AE], with 
E = (ip{0)\H\ip(0)) = (ip(t)\H\ip(t)) = -5.06J. Since 
our systems are small there is generally no meaning to the 
limit AE — > 0, because small enough windows may fail 
to contain even a single eigenstate. Instead, one should 
show that the microcanonical predictions are robust with 
respect to the choice of the width of the energy window. In 
Fig. [6] we demonstrate this robustness in a neighborhood of 
AE = 0.1J, a value that seems to be an appropriate choice 
given the data presented in the inset of the same figure. 
There we show the dependence on AE of the predictions 
for n(k. x = 0) given by the "left-averaged" and the "right- 
averaged" microcanonical ensembles, by which we mean that 
the microcanonical windows are chosen as [Eq ~ AE, Eq] 
and [E ,E + AE], respectively. We see that for AE < 0.1 J, 
the two microcanonical predictions are almost independent 
of the value of AE. The main panel in Fig. [6] shows that the 
microcanonical values of n(k x ) for AE = 0.05J and for 
AE = 0.1J are indistinguishable. 
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FIG. 6: Microcanonical ensemble. Microcanonical momentum 
distribution function for two different values of AE. Inset: Mi- 
crocanonical predictions for n(k x — 0) calculated using the left 
({Eg — AE, Eg]) and the right ([Eo, Eo + AE]) averages as functions 
of AE. 



As the principal observables of interest we chose the 
marginal momentum distribution along the horizontal axis 
n(k x ) = n(k x ,k y ) and in particular its central com- 
ponent n(k x = 0), quantities readily measurable in actual 
experiments with ultracold quantum gases |19[] . Here the 
full two-dimensional momentum distribution is n(k x , k y ) — 

1 / i2 E i ,i e " i27rk(r,_rj)/i <^ t ^>' where L = L x = L y = 5 
are the linear sizes of the lattice. The position r s ; = (i x d , i y d) 
involves the lattice constant d. 

2. The microcanonical ensemble in a small system. 

To compute the microcanonical ensemble predictions, 
we have averaged over all eigenstates whose energies lie 



3. Eigenstate thermalization and the width of the energy 
distribution. 

The eigenstate thermalization alone is not sufficient to en- 
sure an agreement between the predictions of the diagonal 
and thermal ensembles. As discussed in Ref. 11311 . it is also 
necessary that both distributions be sufficiently narrow. More 
specifically, one must require for both ensembles 

(AE) 2 \A"(E)/A(E)\ 1, (8) 

where AE is the width of the energy distribution in the en- 
semble, and A(E) is the dependence of the expectation value 
of the observable A aa — ( , 3> a | J 4|5 , a ) on the energy E a of the 
Hamiltonian-operator eigenstate \^ a )- Note that because of 
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eigenstate thermalization, A(E) is a smooth function of en- 
ergy. For the thermodynamical ensembles the condition (0 is 
always satisfied in the thermodynamic limit. We now show 
that it is also satisfied for the diagonal ensemble in the ther- 
modynamic limit. 

If one considers an observable a that is the intensive coun- 
terpart of A, all conclusions obtained for a can be extended 
to the original observable A via trivial rescaling. For exam- 
ple, for our principal observable of interest, n(k x ), the corre- 
sponding intensive variable is the momentum density £(p x ) 
normalized as jdp x ^(p x ) = 1. Notice that in this case 
f(p*) = n{k x )L x d/(2irN). 

For o, the condition in (0 reduces to 



(A £ )V( e )/a(e)|«l, 



(9) 



where e = E/N. For sufficiently large systems the depen- 
dence of a on e is independent of the system size. Hence, in 
order to justify the validity of (O it is sufficient to prove that 
the width of the distribution of the energy per particle in the 
diagonal ensemble converges to zero for large linear sizes L 
of the system: 



Ae 



0. 



(10) 



Suppose that initially our system is prepared in an eigen- 
state | Wo) of a Hamiltonian Hq and that at time t = the 
Hamiltonian is suddenly changed to H: 

H -> H = H + W, 

where W is the difference between the new and the old Hamil- 
tonians. Within this scenario, the energy width 



AE = 



\ 



a \ a / 



of the diagonal ensemble becomes equal to the variance of the 
new energy in the state |Wq): 



AE = AH = ^(tf |# 2 |*o> - (Wol^l^o) 2 ■ 

It is now straightforward to show that the variance of H equals 
the variance of W: 

AH = AW. 

In order to deduce how AW scales in the thermodynamic 
limit, we assume that W is a sum of local operators w(j) over 
some region of the lattice a (a single point, a straight line, the 
whole lattice, etc.): 

Here w(j) is a polynomial of creation and annihilation oper- 
ators localized at the points j + Aj, where \Aj\ is limited 



from the above by a finite number that does not scale with the 
system size. The mean square of W can be written as 

(*o|^ 2 |*o> = (*o|^|*o> 2 

+ K*oRji)«5(j2)|*o) 
ji,hecr 

- (*oNji)I*o)(*oNj2)|*o)]. (ID 

In the absence of long-range correlations the expression in 
brackets tends to zero for large distances between j\ and fa. 
Therefore, the whole second term on the right-hand-side of 
(fTTT > scales as L d °, where d a is the dimensionality of the sub- 
lattice a and L is the linear size of the lattice. The variance of 
W scales the same way: 



(AWf 



L a 



Retracing our steps, we arrive at the conclusion that the energy 
width Ae indeed tends to zero in the thermodynamic limit: 



a L — >o 

Ae oc 



1 



where dh > d a is the dimensionality of the whole lattice. 

Note that for the two-dimensional lattice considered in this 
paper the role of W is played by the hopping energy of the 
"door" link. An analysis similar to the one above shows that 
increasing the number of "door" links will lead to an increase 
in the width Ae, proportional to the square root of the number 
of "door" links. This is why in our 2D experiment, we have 
chosen the position of the missing sites to be the one in Fig. [5] 
so that only a single link is opened during the time evolution. 
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